STREHL RATIO WITH ANISOPLANATISM

The Strehl ratio (SR) for a circular aperture [#7##1###] from the Huygens--Fresnel approximation is

#math45#

SR = #tex2html_wrap_indisplay859##tex2html_wrap_indisplay860#d#tex2html_wrap_indisplay861#α K(α) exp #tex2html_wrap_indisplay862##tex2html_wrap_indisplay863##tex2html_wrap_indisplay864#.

(1)

#math46##tex2html_wrap_inline866#(#tex2html_wrap_inline867#α) is the structure function, and #math47#K(α) is a factor times the optical transfer function given by

#math48#

K(α) = #tex2html_wrap_indisplay871##tex2html_wrap_indisplay872#cos^-1(α)-α#tex2html_wrap_indisplay873#1-α²#tex2html_wrap_indisplay874##tex2html_wrap_indisplay875# U(1 - α),

(2)

where #math49#U#tex2html_wrap_inline877#x#tex2html_wrap_inline878# is the unit step function defined as

#math50#

U(x)	=	1 for x≥0 ,
U(x)	=	0 for x ;SPMlt; 0 .	(3)

To find the Strehl ratio, one must first determine the structure function. It was found by Fried[#4##1###] for angular anisoplanatism. If the source is collimated and a general displacement is introduced, his expression for a wave propagating from ground to space becomes

#math51#

#tex2html_wrap_indisplay888#(αD)	=	2(2.91) k₀²#tex2html_wrap_indisplay893#dz C_n²(z)#tex2html_wrap_indisplay896#(αD)^5/3+d^5/3(z)
		#tex2html_wrap_indisplay898#-#tex2html_wrap_indisplay899#12 #tex2html_wrap_indisplay900##tex2html_wrap_indisplay901#αD+#tex2html_wrap_indisplay902#d(z) #tex2html_wrap_indisplay903#-#tex2html_wrap_inline904#12#tex2html_wrap_indisplay905# #tex2html_wrap_indisplay906#αD-#tex2html_wrap_indisplay907#d(z) #tex2html_wrap_indisplay908##tex2html_wrap_indisplay909#,	(4)

where #math52#C_n²(z) is the turbulence strength as a function of altitude; #math53#k₀ = 2π/λ, where λ is the wavelength of operation; D is the aperture diameter; and #math54##tex2html_wrap_inline917#d(z) is the vector displacement of the two paths.

The sums of the terms in brackets almost cancel, thus causing difficulties if one tries to evaluate this integral numerically. The terms in the absolute-value sign are equal to

#math55#

#tex2html_wrap_indisplay920# #tex2html_wrap_indisplay921#αD±#tex2html_wrap_indisplay922#d(z) #tex2html_wrap_indisplay923# = #tex2html_wrap_indisplay924##tex2html_wrap_indisplay925#αD#tex2html_wrap_indisplay926#±2αD d (z)cos#tex2html_wrap_indisplay927##tex2html_wrap_indisplay928##tex2html_wrap_indisplay929#+d²(z)#tex2html_wrap_indisplay930#,

(5)

where is the angle between #math56##tex2html_wrap_inline932#α and #math57##tex2html_wrap_inline934#d(z) . This expression can be simplified and the numerical difficulties can be eliminated by using Gegenbauer polynomials.[#8##1###] Their generating function is

#math58#

#tex2html_wrap_indisplay937#1-2ax+a²#tex2html_wrap_indisplay938# = #tex2html_wrap_indisplay939#C_p^λ(x) a^p.

(6)

These functions are sometimes referred to as ultraspherical functions because they are a generalization of the Legendre polynomials P_n(t) , whose generating function is

#math59#

#tex2html_wrap_indisplay945#1-2ax+a²#tex2html_wrap_indisplay946# = #tex2html_wrap_indisplay947#P_p(x) a^p.

(7)

The Gegenbauer polynomials with the cosine of a variable as the argument are given in Eq. (8.934 #2) of Ref. <#93#>8<#93#> and can be rewritten as

#math60#

C_p^λ#tex2html_wrap_indisplay952#cos#tex2html_wrap_indisplay953##tex2html_wrap_indisplay954##tex2html_wrap_indisplay955##tex2html_wrap_indisplay956# = #tex2html_wrap_indisplay957##tex2html_wrap_indisplay958#,

(8)

where - #math61#Γ#tex2html_wrap_inline960#x#tex2html_wrap_inline961# is the gamma function. A particular Gegenbauer polynomial that is required is

#math62#

C₂^-5/6#tex2html_wrap_indisplay966#cos(#tex2html_wrap_indisplay967#)#tex2html_wrap_indisplay968# = #tex2html_wrap_inline969#56#tex2html_wrap_indisplay970#1-#tex2html_wrap_inline971#13cos²#tex2html_wrap_indisplay972##tex2html_wrap_indisplay973##tex2html_wrap_indisplay974##tex2html_wrap_indisplay975#.

(9)

For #math63#αD ;SPMgt; d (z) , the terms in the structure function can be expanded in Gegenbauer polynomials. The zeroth- and all odd-order terms cancel. When the summation index is changed by the substitution #math64#p→2p the result is

#math65#

#tex2html_wrap_indisplay980#(αD) = 2(2.91) k₀²#tex2html_wrap_indisplay983#dz C_n²(z)#tex2html_wrap_indisplay986##tex2html_wrap_indisplay987##tex2html_wrap_indisplay988#.

(10)

It is this canceling of the first two terms of the power series that would cause numerical difficulties. Define a distance moment as

#math66#

d_m≡2.91 k₀²#tex2html_wrap_indisplay993#dz C_n²(z) d^m(z)

(11)

and a phase variance as

#math67#

σ_{#tex2html_wrap_indisplay1000#}² = d_5/3.

(12)

Unlike the calculation for Strehl ratio for uncorrected turbulence and for corrected turbulence with tilt jitter, an exact analytical solution cannot be found for anisoplanatism. Fortunately, for adaptive-optics systems, the Strehl ratio should be fairly high by design, which requires the structure function to be small. This assumption allows one to retain only the first term of the Gegenbauer expansion to give

#math68#

#tex2html_wrap_indisplay1003#(αD) = 2σ_{#tex2html_wrap_indisplay1006#}² - 2x,

(13)

where

#math69#

x = d₂#tex2html_wrap_indisplay1009#1-#tex2html_wrap_inline1010#13cos²#tex2html_wrap_indisplay1011##tex2html_wrap_indisplay1012##tex2html_wrap_indisplay1013##tex2html_wrap_indisplay1014##tex2html_wrap_indisplay1015#56(αD)^-1/3.

(14)

We justify this single-term approximation below by showing that it produces a result close to the exact result.
...
The Strehl ratio with the six term approximation is

#math70#

SR #tex2html_wrap_indisplay1018# #tex2html_wrap_indisplay1019##tex2html_wrap_indisplay1020#d#tex2html_wrap_indisplay1021#α K(α) #tex2html_wrap_indisplay1022##tex2html_wrap_indisplay1023##tex2html_wrap_indisplay1024#.

(15)

If just the first term in the last parenthetical expression is retained, the result is equivalent to the extended Maréchal approximation. It is shown below that the six-term approximation is best for aperture sizes normally encountered. The angle integral for the nth term, after use of the binomial theorem, is proportional to

#math71#

Φ(n) = #tex2html_wrap_indisplay1028##tex2html_wrap_indisplay1029#d#tex2html_wrap_indisplay1030# #tex2html_wrap_indisplay1031#1-#tex2html_wrap_indisplay1032#13cos²#tex2html_wrap_indisplay1033##tex2html_wrap_indisplay1034##tex2html_wrap_indisplay1035##tex2html_wrap_indisplay1036# = #tex2html_wrap_indisplay1037##tex2html_wrap_indisplay1038##tex2html_wrap_indisplay1039##tex2html_wrap_indisplay1040##tex2html_wrap_indisplay1041# 3^-m#tex2html_wrap_indisplay1042#d#tex2html_wrap_indisplay1043# cos^2m#tex2html_wrap_indisplay1044##tex2html_wrap_indisplay1045##tex2html_wrap_indisplay1046#,

(16)

where

#math72#

#tex2html_wrap_indisplay1049##tex2html_wrap_indisplay1050##tex2html_wrap_indisplay1051# = #tex2html_wrap_indisplay1052#.

(17)

Equation (4.641 # 4) in Gradshteyn and Ryzhik[#8##1###] is

#math73#

#tex2html_wrap_indisplay1055#d#tex2html_wrap_indisplay1056# cos^2m#tex2html_wrap_indisplay1057##tex2html_wrap_indisplay1058##tex2html_wrap_indisplay1059# = #tex2html_wrap_indisplay1060#,

(18)

where

#math74#

(2m - 1)!!	=	(2m - 1)(2m - 3)…(3)(1),	(19)
(2m)!!	=	(2m)(2m - 2)…(4)(2).	(20)

With these relations, the angle integral is equal to

#math75#

Φ(n) = 1 - #tex2html_wrap_indisplay1070##tex2html_wrap_indisplay1071##tex2html_wrap_indisplay1072##tex2html_wrap_indisplay1073# 3^-m#tex2html_wrap_indisplay1074#.

(21)

The values of interest to us are #math76#Φ(0) = 1, #math77#Φ(1) = 0.8333, #math78#Φ(2) = 0.7083, #math79#Φ(3) = 0.6134, #math80#Φ(4) = 0.5404, and #math81#Φ(5) = 0.4836. The aperture integration for the nth term is proportional to

#math82#

Y(n) = #tex2html_wrap_indisplay1084#dα α^1-n/3K(α).

(22)

This is a generalization of the integral evaluated by Tatarski in Sec. 55, Eq. (22) of Ref. <#217#>9<#217#>. Its value is

#math83#

Y#tex2html_wrap_indisplay1087#n#tex2html_wrap_indisplay1088# = #tex2html_wrap_indisplay1089# Γ #tex2html_wrap_indisplay1090##tex2html_wrap_indisplay1091##tex2html_wrap_indisplay1092# ;SPMnbsp;;SPMnbsp;;SPMnbsp;;SPMnbsp; for n ;SPMlt; 6.

(23)

The values of interest to us are #math84#Y#tex2html_wrap_inline1094# 0#tex2html_wrap_inline1095# = 1, #math85#Y#tex2html_wrap_inline1097#1#tex2html_wrap_inline1098# = 1.402, #math86#Y#tex2html_wrap_inline1100#2#tex2html_wrap_inline1101# = 2.087, #math87#Y#tex2html_wrap_inline1103#3#tex2html_wrap_inline1104# = 3.396, #math88#Y#tex2html_wrap_inline1106#4#tex2html_wrap_inline1107# = 6.419, and #math89#Y#tex2html_wrap_inline1109#5#tex2html_wrap_inline1110# = 16.94. With these values for the integral, the Strehl ratio approximation is

#math90#

SR #tex2html_wrap_indisplay1113# (1+0.9736 E+0.5133 E²+0.2009 E³+0.0697 E⁴+0.02744 E⁵)exp(-σ_{#tex2html_wrap_indisplay1116#}²),

(24)

where

#math91#

E = #tex2html_wrap_indisplay1119#.

(25)

...
There is an error made in using this approximation for the central part of the aperture that increases with each term in the approximation. One has to determine whether this error is less than or greater than the increased accuracy achieved over the remainder of the aperture by using additional series terms. To resolve these uncertainties, I compared the Strehl ratio, using various numbers of terms, with exact calculations.

I calculated the Strehl ratio numerically for the case in which the displacement does not vary with propagation distance. In Fig.~#f1#236> are plotted the exact Strehl ratio versus displacement for the Hufnagel--Valley 21 (HV-21) model of turbulence[#10##1###,#11##1###,#12##1###] and the Strehl ratio from relation (24) for #math92#D/r_o = 1, with only the unity term in parenthesis (extended Marechal approximation) and with different numbers of terms in the parenthesis.
...